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ABSTRACT 

Several recently discovered peculiar type la supernovae seem to demand an altogether new formation 
theory that might help explain the puzzling dissimilarities between them and the standard type la 
supernovae. The most striking aspect of the observational analysis is the necessity of invoking super- 
Chandrasekhar white dwarfs having masses ~ 2.1 — 2.8Mq, Mq being the mass of Sun, as their most 
probable progenitors. Strongly magnetized white dwarfs having super- Chandrasekhar masses were al- 
ready established to be potential candidates for the progenitors of peculiar type la supernovae. Owing 
to the Landau quantization of the underlying electron degenerate gas, theoretical results yielded the 
observationally inferred mass range. Here we sketch a possible evolutionary scenario by which super- 
Chandrasekhar white dwarfs could be formed by accretion on to a commonly observed magnetized 
white dwarf, invoking the phenomenon of flux freezing. This opens the multiple possible evolutions 
ending in supernova explosions of super-Chandrasekhar white dwarfs having masses within the range 
stated above. We point out that our proposal has observational support, like, the recent discovery of 
a large number of magnetized white dwarfs by SDSS. 

Subject headings: accretion, accretion disks — equation of state — novae, cataclysmic variables — 
stars: magnetic field — supernovae: general — white dwarfs 



1. INTRODUCTION 

Identifying the progenitors of type la supernovae is 
an extremely important and ongoing research issue 
(iHowell 20111 ). These supernovae are considered to be 
standard candles for cosmic distance measurements and 
hence are pr obes for studying the expansion history of 
the universe (jPerlmutter et al. 19 99). According to the 
general consensus, a carbon-oxygen white dwarf in a bi- 
nary system accretes matter from a companion star and 
gradually approaches the Cha ndrasekhar mass limit of 
1.44Mq ijChandrasekhar 1935f) . M© being the mass of 
Sun. This further triggers a thermonuclear instability 
in the white dwarf, resulting in a violent and extremely 
energetic explosion, which we observe as a type la super- 
nova. However, far less is known about the nature of the 
companion star (jParthasarathv et al. 20071 ). 

Adding to the puzzle further is the recent discovery of 
several peculiar type la supe rnovae — SN 20 06gz, SN 
2007if, SN 2009dc, SN 2003fg (|Scalzo et al. 20101 . These 
supernovae are distinctly over luminous, powered by a 
higher than usual production of nickel, and have very low 
ejecta velocity com pared to their standard counterparts 
(jHowell et al. 2006r). Thev also violate the luminosity- 
stretch relation (jPhillips 19931 : iGoldhaber et al. 20011 ) 
which prohibits them from being categorized as 
standard candles. Interestingly, however, these 
anomalies seem to be resolved if instead of the 
standard theory one invokes super- Chandrasekhar- 
mass white d warfs, with masses 2.1 — 2.8M(t,, as their 

progenitors (fScalzo et al. 20101 : Howell ct al. 2006; 
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[Silverman et al. 20111 : iTaubenberger et al. 20l"ll ). 

The leading question now is the origin and stability 
of these super-Chandrasekhar white dwarfs. There 
have been a few simulations of accreting binary white 
dwarfs, whi ch include differential rotation and other 
parameters ()Hachisu et al. 2012) . that try to explain 
the required mass range stated above. In addition, 
models have been proposed involving mergers of two 
white dwarfs, namely t he double degener ate scenario 
(|Iben fc Tutukov 1984 IScalzo et al. 20101) . and core 
degenerate scenario ( Kashi fc Soker 20111 ). However, a 
detailed theoretical understanding seems to be lacking 
in them. 

On a fundamentally different ground. Das & 
Mukhopadhyay (2012a, 2012b) have shown that strongly 
magnetized white dwarfs, having central fields ~ 10^^ — 
10^^ G, are capable of having super-Chandrasekhar 
masses in the range 2 — 2.6Mq. The basic idea employed 
(Das & Mukhopadhyay 2012a, 2012b, hereafter DM1, 
DM2 respectively) is that, for magnetic fields (B) greater 
than a critical field Be = 4: AM x 10^^ G, the effect of Lan- 
dau quantization on the underlying electron degenerate 
matter becomes significant and as a result the densit y 
of states for the electrons changes ()Lai fc Shapiro 19^11) . 
This in turn modifies the corresponding equation of state 
(EoS) and hence the mass-radius relation of these white 
dwarfs. Now the additional question is to form white 
dwarfs with such strong magnetic fields. 

In this Letter, we sketch an evolutionary scenario by 
which highly magnetized white dwarfs could be formed. 
We concentrate on the single degenerate progenitor sce- 
nario. Starting with some reasonable assumptions about 
the mass and the magnetic field of a normally observed 
magnetic white dwarf, we demonstrate that by the mech- 
anism of flux freezing and Landau quantization a possi- 
ble evolutionary link between strongly magnetized super- 
massive white dwarfs and the observed peculiar type la 
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supernovae could be established. This scenario is com- 
patible with the discovery of several (isolated) magne- 
tized white dwarfs by SPSS having; surface fields 10^ — 
10^ G (|Schmidt et ah 20031 IVmlandingham et al. 20051) 
and the fact that 25% of the observed cataclysmic vari- 
ables (CVs) have surface field strengths as high as 
10^ - 10* G (jWickramasinghe fc Ferrario 20001) . 

This Letter is organized as follows. In §2 we describe 
how a magnetized, sub-Chandrasekhar, accreting white 
dwarf evolves into a super- Chandrasekhar white dwarf. 
We subsequently discuss the steps that lead to the forma- 
tion of multiple evolutionary tracks in the mass-radius 
plane, which in turn help in explaining the observed 
super-Chandrasekhar (progenitor) mass range of the pe- 
culiar type la supernovae. In §3 we study the timescale 
of evolution of these magnetized accreting white dwarfs. 
Finally we conclude in §4 with a discussion about the 
observational evidence supporting our proposed evolu- 
tionary scenario. 

2. EVOLUTIONARY PATH FROM A 
SUB-CHANDRASEKHAR TO A 
SUPER-CHANDRASEKHAR WHITE DWARF 
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Fig. 1. — Mass- Radius relations — solid line represents Chan- 
drasekhar's non-magnetic result; dotted, dashed and dot-dashed 
lines represent the evolutionary tracks of accreting magnetized 
white dwarfs, ending in super-Chandrasekhar type la supernovae 
with progenitor masses 2.01Mq, 2.33Mq and 2.58M0 respectively. 
The solid lines marked with squares represent from top to bot- 
tom the intermediate mass-radius relations corresponding to 50124- 
level, 200-level and 1-level systems of Landau quantization respec- 
tively. 

We consider a white dwarf having a central field sev- 
eral orders of magnitude higher than that of the sur- 
face, with a typical mass and radius to be Mq ~ O.2Af0 
and R ^ 15000 km respectively. Note, as justified in 
§IV.C of DM1, that the magnetic field might be approx- 
imately constant in a range of radii around the center 
of the white dwarf, which we define as internal mag- 
netic field. This is quite a plausible assumption, since 
the original star collapsing to form a white dwarf might 
have a very large interior magnetic field compared to 
that obse rved on its surface, as h as been the case for 
the Sun (i Gough fc Mclntvre 199 ^. Thus, it is quite 
likely that the values of the central and surface mag- 
netic fluxes (in the initial star and hence in the white 



dwarf) would be significantly different, developing a core- 
envelope boundary. Hence, we start with an internal 
magnetic field Bint ~ 5 x 10^^ G along with a surface 
field Bs at least three orders of magnitude lower than 
Bint- Such a white dwarf still has B < Be and it lies 
on Chandrasekhar 's (non- magnetic) mass-radius relation 
having EoS similar to that of a nonmagnetic case. Now 
as the white dwarf accretes matter, there is an amplifi- 
cation of its central magnetic field as a consequence of 
the increase in central density (via flux freezing theo- 
rem) due to the contraction in size of the white dwarf 
(|Cumming 2002( ). This leads the central magnetic field 
to eventually exceed Be- Hence, EoS of the underlying 
electron degenerate gas and the mass-radius relation of 
the underlying white dwarf modify, as shown in DM1 
and DM2. As a result, the white dwarf transits from the 
mass-radius curve for the nonmagnetic EoS to that for a 
magnetic EoS. As the outward modified pressure counter- 
acts the inward gravitational force, a quasi-equilibrium 
state attains, which is determined by the degree of Lan- 
dau quantization of the system. Larger the magnetic 
field, smaller is the number of Landau levels occupied 
by the electrons, leading to more massive white dwarfs 
(DM1, DM2). Hence, as accretion continues, the mass- 
radius curve of the white dwarf deviates more and more 
from that of Chandrasekhar's. In Figurefl] we show three 
mass-radius relations (solid curves marked with squares) 
following DM1, for magnetized white dwarfs with three 
different values of Bint — 1-7 x 10" G, 4.4 x lO^^ G 
and 8.8 x 10^^ G, for the top to bottom curves respec- 
tively, through which the evolving white dwarf, with ini- 
tial mass, radius and magnetic fields mentioned above, 
passes at different intermediate stages. 

Following Bandyopadhyay et al. (1997), we adopt a 
profile for magnetic field in our study assuming the sur- 
face and central magnetic fluxes of the initial and final 
white dwarfs are individually conserved as the evolu- 
tion progresses. For example, for the white dwarf with 
B,nt = 8.8 X 10^^ G and i? = 69.5 km (which is the 
maximum possible radius for the corresponding mass- 
radius relation), if we consider Bint to be confined upto 
Rint — R/IO (corres ponding to a profile with par ameters 
/3 = 0.01,7 = 3.6 (|Bandvopadhvav et al. 1997t )). then 
the above assumption holds true with surface and central 
fiuxes respectively 1.46 x 10^'' G cm^ and 4.25 x 10^^ G 
cm^. Note, however, that the extent of Bint adjusts ac- 
cordingly as the white dwarf evolves and hence Rint need 
not be the same for the initial and final white dwarfs. 
Figure [1] shows how the an initial mass of the white 
dwarf could possibly evolve to 3 different final limiting 
masses, namely 2.OIM0, 2.33M0 and 2.58il/0. Once 
the corresponding final mass is reached, further accretion 
induces a runaway thermonuclear reaction, causing the 
white dwarf to explode at that mass, leading to a super- 
Chandrasekhar supernova. We note that, although we 
show only three tracks in Figure [U any number of tracks 
are possible resulting in different final limiting masses 
< 2.58M0. In the following subsection we explain in 
detail how we obtain these multiple tracks. 

2.1. Multiple evolutionary tracks 

When the electrons occupy the ground Landau level, 
pressure increases monotonically as a function of den- 
sity. This corresponds to a mass-radius relation where 
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initially the mass increases with increasing radius (let us 
call it the first branch of the mass-radius relation). If 
the magnetic field is sufhciently high, then the radius 
becomes nearly independent of mass at higher central 
densities in this branch. From DM1, we note that when- 
ever there is a transition from a lower Landau level to the 
next higher level, a kink appears in EoS, followed by a 
plateau — which is a small region where the pressure be- 
comes nearly independent of the density (see Figure 1 of 
DM1). For a given magnetic field, the central density of 
the white dwarf corresponding to the kink in EoS arising 
at the transition from the ground to first Landau levels 
corresponds to a white dwarf having the maximum possi- 
ble radius. Note that the systems having more than one 
occupied Landau level exhibit multiple branches in their 
mass-radius relations including the 'first branch' (for de- 
tails see DM1). For example, the mass-radius relation 
of white dwarfs having central densities corresponding 
to the plateau of EoS shows a turning point and subse- 
quent decrease of mass with decreasing radius (see Figure 
2 of DM1). 

Now, let us consider a white dwarf lying on an interme- 
diate mass-radius relation in its early stage of evolution 
(e.g. the top most solid curve marked with squares in 
Figure [T]), such that it has the maximum possible ra- 
dius. As it evolves via accretion, its magnetic field and 
mass increases with decreasing radius, such that when 
the white dwarf transits to the next mass-radius relation, 
it still has the corresponding maximum radius satisfying 
the flux-freezing condition. For a sufficiently high field, 
the white dwarf with maximum radius also corresponds 
to the maximum possible mass. Thus, in Figure 1, we 
show that the evolution of such a white dwarf gives rise to 
the track corresponding to the limiting mass of 2.58M0, 
which corresponds to the new mass limit for white dwarfs 
proposed by Das & Mukhopadhyay (2013). Hence, this 
track demarcates a zone in the mass-radius plane, beyond 
the right-hand side of which no further track is possible. 

From the bottom two mass-radius relations (represent- 
ing further intermediate masses and radii of the evolving 
white dwarf) in Figure 1, we notice that the evolving 
white dwarf passing through them may lie anywhere in a 
zone where the radius is independent of mass satisfying 
the flux freezing condition. This zone is a part of the 'first 
branch' of the corresponding mass-radius relation. This 
happens at a sufficiently high field > 1.7 x 10"'^"' G (which 
corresponds to a 5000-Landau- level system), when a part 
of EoS can be described by a polytropic relation of the 
form 

P = i^/, when r = 2, (1) 

where P and p are respectively the pressure and density 
of the electron gas, F is the polytropic index and K is 
a constant depending on the magnetic field. Now on 
solving the magnetostatic equilibrium condition by Lane- 
Eniden formalism, one obtains the following scaling laws 
(DM1) for radius (i?) and mass {M): 

R oc (2) 

and 

M oc p(3r-4)/2^ (3) 

where pc is the central density of the white dwarf. Thus 
when r = 2, the radius becomes independent of the 
central density, while the mass becomes proportional to 
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the central density. This indicates that the same mag- 
netic field leads to a white dwarf having different possible 
masses with the same radius. Hence once the above poly- 
tropic relation starts holding true, the fiux freezing con- 
dition opens the multiple possible evolutionary tracks. 

3. TIMESCALE OF MASS EVOLUTION 

Here we explore the typical timescale of evolution of 
the accreting magnetic white dwarf systems with as- 
sumed parameters quite close to the class of magnetic 
CVs called the Intermediate Polars (IPs). 

3.1. Constant mass accretion rate 
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Fig. 2. — Mass evolution of accreting magnetized white dwarfs 
— dot-dashed, dashed and dotted lines represent constant accre- 
tion rates of respectively 2.38 X IQ-^MQyr-l, 2.15 X lO-^M0yr-l 
and 1.83 X 10~^MQyr~^ for supernova explosions occurring at the 
same time, for progenitor masses 2.58Mq, 2.33M0 and 2.01^0 
respectively. 

Figure [2] shows the evolution of an initial white dwarf 
of mass O.2M0 and radius 15000 km into the 3 super- 
Chandrasekhar white dwarfs with M > 2Mq corre- 
sponding to Figure [TJ Here we assume a constant mass 
accretion rate (M), represented by the constant slopes 
of the 3 lines in Figure [2l for each case, such that they 
explode at the same time. The values of M (typical of 
IPs) are 2.38 x lO-^Mgyr-i, 2.15 x IQ-^Moyr-i and 
1.83 X IQ-^Moyr-i for the exploding masses of 2.58M0, 
2.33Mq and 2.OIM0 respectively. We see that the sys- 
tem with a higher M accumulating more mass in the 
same time (roughly a billion year) gives rise to a more 
massive supernova explosion. 

Figure 3(a) shows that how the exploding mass of 
2.33M0 is attained by 3 different Ms — 3x lO-^Moyr"!, 
2.5 X lO-^Mgyr^i and 2 x lO^^Afoyr^i ~ resulting in 
the explosions occurring at different times. 

3.2. Varying mass accretion rate 

The most favored models for type la supernova 
progenitors invoke a very high AI > lO^^AfQyr"^, 
such that the accreted hydrogen and helium can burn 
steadily ijCumming 2002[ ). These values of M are ob- 
served in supersoft X-ray sources and symbiotic binaries 
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Fig. 3. — Mass evolution of accreting magnetized white dwarf 
with (a) constant accretion rate — sohd, dashed and dotted 
Unes represent 3 X IQ-^MQyr-l, 2.5 X IQ-^MQjr-'^ and 2 X 
lO~^M0yr~^ respectively, (b) varying accretion rate — solid and 
dashed lines correspond to « = 2 and a = I respectively. 

()Cuinniing 20021: iParthasarathv et al. 20071 ). At lower 
Ms, hydrogen burning is unstable and occurs in flashes, 
while at higher rates an ex tended envelope is formed 
(jvan den Heuvel et al. 1992[ ). 

We now plan to invoke a range of Ms, low to high, in 
a single model evolution of the white dwarf and hence 
consider 

where a is a parameter that determines the functional 
dependence of M on the instantaneous mass {Al > Mq) 
of the white dwarf. 

In Figure 3(b) we show two cases corresponding to 
a — 1 and a — 2, having the same final exploding mass, 
namely 2.33M0. In the case with a — 1, M varies from 
lO~^M0yr-i to 1.3 x lO'^Mgyr'^ The evolution is 
slow in this case and the supernova explosion occurs in 
~ 5 X 10* years. In the case with a = 2, M varies from 
IQ-^Msyr-i to 1.7 x IQ-'^Mgyr-^ The evolution here 
is more than twice as fast as the previous case and the su- 
pernova explosion occurs in ~ 2 x 10* years. Hence, a > 2 
is supposed to validate the underlying white dwarfs as 
progenitors of the observed peculiar type la supernovae. 

4. DISCUSSION AND CONCLUSIONS 

Based on the theoretical work on strongly magnetized 
white dwarfs having super-Chandrasekhar masses (DM1, 
DM2), we have discussed possible evolutionary paths 
that lead from magnetized accreting white dwarf bina- 
ries to the peculiar type la supernova explosions. Mul- 
tiple evolutionary tracks are possible in the mass-radius 
plane, which covers a mass range ^ 2 — 2.6Mq, that lies 
within the observational limits. We have considered both 
constant as well as varying Mq, in order to estimate the 
timescale of occurrence of the supernova events. 

Highly magnetized white dwarfs have been discov- 
ered b y SDSS having surface fields in the range 10^ — 
10^ G ( Schmi dt et al. 2 003: Va nlandingham et al. 2001 
IWickramasinghe &: Ferrario 2000 t). Note that their cen- 
tral fields are expected to be a few orders of magnitude 



higher. The properties that we have assumed in our com- 
putation are similar to those seen in IPs. They have 
surface magnetic fields ranging 10^ — 10^ G and they 
are high accre tors with M « (0.2 — 4) x lO^^Moyr"^ 
(jWarner 19951 ). Gumming (2002) studied the evolution 
of magnetic fields in accreting white dwarfs and found 
that surface magnetic fields are reduced significantly for 
M > Mc w (1 - 5) x lQ-i°MQyr-\ due to the advec- 
tion of the field into the interior of the white dwarf by 
the accretion flow. This presumably increases the cen- 
tral magnetic field (jGumming 20021 ). Further, certain 
dwarf novae (a subclass of CVs) like RU Peg, are ex- 
pected to show short and long term variations in Ms, 
ranging from as l ow as (1 — 2) x 10~^M(^yr~ ^ to as high 
as lO~'^M0yr"i (jParthasarathv et al. 20071) . Hence, it 
is quite conceivable that IPs with high internal field and 
low surface field with an increasing M can lead to super- 
Ghandrasekhar white dwarfs causing the peculiar super- 
nova la explosions. Moreover, for a = 2 according to our 
chosen model, the desired M for a type la supernova ex- 
plosion to occur is attained starting from an initial, lower 
rate, typically observed in IPs. Note that peculiar type la 
supernovae arise, according to our argument, from highly 
magnetized GVs, which are presumably 10% of all mag- 
netized GVs. Hence about 2.5% of all GVs should lead 
to super-Ghandrasekhar objects, which is consistent with 
the observed rate (1 - 2%) of peculiar type la supernovae 
(jScalzo et al. 2012D . 

Apart from the advection of the magnetic field into the 
interior caused by accretion, there are other ways of gen- 
erating such strong fields.'^ Interestingly, recent observa- 
tions of magnetic Ap and Bp stars indicate that they have 
a magnetic flux in the range 10^^ — 10^^ G cm^, which is 
very similar to that observed in ma gnetized white dwarfs 
(jFerrario fc Wickramasinghe 20051 ) . This strongly points 
towards the fact that the magnetic fields of highly mag- 
netized white dwarfs are fos sil remnants from their main- 
sequence progenitor stars ()Tout et al. 20041 ). Simula- 
tions of magnetic field evolution in Ap stars strongly sup- 
port the idea of this field itself being a fossil remnant from 
the interstellar medium, generated during the process of 
star formation (Braithwaitc & Spruit 2004). Results of 
these simulations to hold true, the new born star is re- 
quired to have an initial field mainly confined to the core, 
which can be justi fied from the flux freezing the ory dur- 
ing its formation ( Braithwaite &: Nordlund 20061 ). More- 
over, the contraction of a typical interstellar cloud of ra- 
dius ~ 0.1 pc, mass ^ Mq and having a frozen-in mag- 
netic field ~ 3 X 10^^ G could give rise to a fie ld ~ 10* 
G in the resulting star (jShapiro fc Teukolskv 19 83) with 
solar radius. Thus, a star with such a high centrally con- 
centrated field is a plausible progenitor of the strongly 
magnetized super-Ghandrasekhar white dwarfs discussed 
in this Letter. These stars also have convective cores and 
in some special cases a dynamo mechanism might occur 
to generate a high magnetic field therein, leading to a 
(central) magnetic flux as high as ^ 10'^'^ G cm^, sim- 
ilar to that of the magnetized white dwarfs mentioned 
in §2. The incidence of such main-sequence magnetic 
stars might be rare, but so are the super-Ghandrasekhar 

^ However, see the Appendix of DM1, which discusses that even 
a relatively weaker field than that considered here can give rise to 
the super-Chandrasekhar masses. 
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supernova events. 

If the single-degenerate accreting scenario is the correct 
progenitor model for the overluminous (peculiar) type 
la supernovae, then observationally such high magnetic, 
high mass white dwarfs with small radius (also with high 
M) could be one of the plausible progenitors of these su- 
pernovae. However, a detailed investigation needs to be 
carried out in order to understand the physics of explo- 
sions of these white dwarfs and if they indeed produce 
the light curves similar to that of the peculiar type la 
supernovae. 

Now the white dwarfs would be seen as peculiar, super- 
Chandrasekhar objects only during a small fraction of 
their life time, which is ^ 0.05 for a — 2. Hence, 
about ~ 0.12% of the CVs with varying M would be 
likely candidates for peculiar objects. The varying M is 



further justified by astronomical indications, which im- 
ply higher accretion rates at later ti mes of the evolution 
of magnetized mass i ve white dwarfs (I Wang Sz Han 20121 : 
iToonen et al. 2013t iNomoto et al. 2007D . We specu- 
late that these white dwarfs will have X-ray luminosity 
in between commonly observed accreting white dwarfs 
and accreting neutron stars (in the range of 10^^ - 
Iq38 gj.gg g-i). Bright X-ray selected CVs (and other 
objects) might harbor some of these peculiar super- 
Chandrasekhar objects. 



This work is partly supported by ISRO project Grant 
No. ISRO/RES/2/367/10-11. U.D. thanks CSIR, India 
for financial support. The authors would like to thank 
Arnab Rai Choudhuri of IISc for discussion. 



REFERENCES 



Bandyopadhyay, D., Chakrabarty, S., & Pal, S. 1997, Phys. Rev. 

Lett., 79, 2176. 
Braithwaite, J., & Spruit, H. C. 2004, Nature, 431, 819. 
Braithwaite, J., & Nordlund, A. 2006, A&A, 450, 1077. 
Chandrasekhar, S. 1935, MNRAS, 95, 207. 
Chanmugam, G., Ray, A., & Singh, K. P. 1991, ApJ, 375, 600. 
Gumming, A. 2002, MNRAS, 333, 589. 

Das, U., & Mukhopadhyay B. 2012a, Phys. Rev. D, 86, 042001; 
DM1. 

. 2012b, Int. J. Mod. Phys. D, 21, 1242001; DM2. 

Das, U., & Mukhopadhyay, B. 2013, Phys. Rev. Lett., 110, 071102. 

Ferrario, L., & Wickramasinghe, D. T. 2005, MNRAS, 356, 615. 

Goldhaber, G., et al. 2001, ApJ, 558, 359. 

Gough, D. O., & Mclntyre, M. E. 1998, Nature, 394, 755. 

Hachisu, I., Kato, M., Saio, H., & Nomoto, K. 2012, ApJ, 744, 69. 

Hicken, M., et al. 2007, ApJ, 669, L17. 

Howell, D. A. 2011, Nature Gomm., 2, 350. 

Howell, D. A., et al. 2006, Nature, 443, 308. 

Iben, I., & Tutukov, A. V. 1984, ApJS, 54, 335. 

Kashi, A., & Soker, N. 2011, MNRAS, 417, 1466. 

Lai, D., & Shapiro, S. L. 1991, ApJ, 383, 745. 

Nomoto, K., Saio, H., Kato, M., & Hachisu, I. 2007, ApJ, 663, 
1269. 

Parthasarathy, M., Branch, D., Jeffery, D. J., & Baron, E. 2007, 
New Ast. Rev., 51, 524. 



Perlmutter, S., et al. 1999, ApJ, 517, 565. 
Phillips, M. M. 1993, ApJ, 413, L105. 
Scalzo, R. A., et al. 2010, ApJ, 713, 1073. 
Scalzo, R. A., et al. 2012, ApJ, 757, 12. 
Schmidt, G. D., et al. 2003, ApJ, 595, 1101. 

Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs 
and Neutron Stars: The physics of compact objects (John Willey 
& Sons, New York). 

Silverman, J. M., et al. 2011, MNRAS, 410, 585. 

Taubenberger, S., et al. 2011, MNRAS, 412, 2735. 

Toonen, S., Nelemans, G., Bours, M., & Portegies Zwart, S. 
larXiv:1302 .0837 

Tout, G. A., Wickramasinghe, D. T., & Ferrario, L. 2004, MNRAS, 
355, L13. 

van den Heuvel, E. P. J., Bhattacharya, D., Nomoto, K., & 

Rappaport, S. A. 1992, A&A, 262, 97. 
Vanlandingham, K. M., et al. 2005, AJ, 130, 734. 
Wang, B., & Han, Z. 2012, New Ast. Rev., 56, 122. 
Warner, B. 1995, Gataclysmic Variable Stars (Gambridge: 

Gambridge Univ. Press) 
Wickramasinghe, D. T., & Ferrario, L. 2000, Pub. Ast. Soc. Pac, 

112, 873. 

Yamanaka, M., et al. 2009, ApJ, 707, L118. 



